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Abstract. - The fluctuations in the spacing of the tunneling resonances through a quantum 
dot have been studied in the quantum Hall regime. Using the fact that the ground-state of 
the system is described very well by the Laughlin wavefunction, we were able to determine 
accurately, via classical Monte Carlo calculations, the amplitude and distribution of the peak- 
spacing fluctuations. Our results clearly demonstrate a big enhancement of the fluctuations 
as the importance of the electronic correlations increases, namely as the density decreases and 
filling factor becomes smaller. We also find that the distribution of the fluctuations approaches 
a Gaussian with increasing density of random potentials. 



Several experimental studies |l|, |2|, |3j have recently demonstrated that the fluctuations in 
the ground-state energy of a quantum dot, which are manifested in the fluctuations in the 
resonant-tunneling-peak spacings, are much larger than what one expect from models that 
ignore electron correlations. Numerical studies Jl], |], ||, || have indeed revealed an enhancement 
of the ground-state energy fluctuations due to electron-electron interactions. 

In this work we present calculations for an interacting electron system in a regime that can 
be treated almost exactly - the quantum Hall regime. The ground-state wavefunction in this 
regime is faithfully described by the Laughlin wavefunction 0. Consequently, as long as the 
potential fluctuations do not mix in excited states (i.e. when the potential energy is smaller 
than the gap), the peak-spacing fluctuations (PSF) can be evaluated as expectation values 
in the Laughlin state. As such expectation values can be easily calculated using the plasma 
analogy |7), via, e.g., classical Monte Carlo simulations, we are able to obtain the magnitude 
of the PSF, their distribution and their dependence on the range of the potentials and the 
electron number. Indeed we find that the more important the electronic correlations (the 
lower the filling factor), the larger the magnitude of the PSF. In addition, we also find that 
the distribution of PSF is Gaussian, in agreement with experiments jl], ^, || . 

The experimentally measured quantity is the spacings between the resonant-tunneling peaks 
through the quantum dot. At low enough temperatures (smaller than the excitation energies 
of the dot), the peak spacing is determined by the addition spectrum, 



model § Ep = N(N - l)U/2 + £f where U is the char ging energy and ej is the 




Typeset using EURO-TeX 



2 



single-electron spectrum. In this model the peak-spacing fluctuations are determined by the 
fluctuations in the single-electron spectrum, which are described by random-matrix theory, 
in contrast with the experimental observations. This deviation from random-matrix theory 
was attributed to the importance of electronic correlations, which are not captured in the 
constant-interaction model. 

For strong magnetic fields, in particular in the fractional quantum Hall regime, electronic 
correlations become very significant and dominate the underlying physics. The ground-state 
wavefunction of ./V electrons in a quantum dot in a strong magnetic field, with filling factor 
v = 1/m, can be described very well by the Laughlin wavefunction |Q 

( Z1 , z N ) = n ^ - *iT e " E - |z,|2/4 > ( 2 ) 

i<j 

where Zi denotes the complex coordinates of the i-th particle, Zj = Xi +iyi, and all lengths are 
expressed in units of the magnetic length, £jj = ^hc/eH. 

Since the magnetic field quenches the kinetic energy and the interaction energy varies 
smoothly with the number of particles, the fluctuations in the peak spacings stem only from 
fluctuations in the potential energy. If the energy gap due to correlations is large enough 
(compared to the potential), the potential energy in the presence of a random potential V{z), is 
given by £S = J d 2 rp^ /m) (r)V( r ), where p« /m) (r) is the single-particle distribution function, 

P% /m) (r) = j d 2 z 1 ...d 2 z N \¥ 1 ' m) {z 1 ,...,z N )\ 2 5{\ Zl \ 2 -r 2 ) (3) 

The single-particle distribution p { ^, /m) (r) can be calculated using the mapping onto a classical 
plasma model Q. This mapping is based upon rewriting l^j 2 as 



exp 



^\zi\ 2 /2 + m^log(\zi 



(4) 



Consequently, expectation values in the ground-state can be expressed as statistical aver- 
ages for a classical system of N particles in a harmonic confining potential and logarithmic 
interactions. These averages were evaluated using a classical Monte Carlo approach |J. 

In fig. la we plot the single-particle density of 15, 20 and 25 particles in the v = 1/3 
fractional quantum Hall regime. As can be seen from the figure, the bulk density is constant, 
while the edge structure for different particle numbers is identical, but relatively shifted. Since 
the edge of the iV-particle system in the v = 1/m case is determined by the maximal occupied 
angular momentum state, m(N — 1), localized at distance y / 2m(N — 1) from the origin, the 
edge structure of N% and ^-particle systems can be made to collapse by a relative shift 
of y/2m(Ni — 1) — ^/2m(N2 — 1). In fig. lb we replot the same densities as in fig. la, but 
relatively shifted in that manner. Indeed we see that all the distributions collapse onto a single 
curve. This invariance of the edge-structure allows us to deduce immediately the TV-dependence 
of the peak-spacing fluctuations: writing 



A(«) _ 
2 — 



d A r A 2 p£'"»(r) V( r) 



d 2 r 



dN 2 



V(r), 



(5) 



(with &2pN = Pjv+i + Pn-i ~~ ^Pn), and using the abovementioned fact that 



/^ /m) (r) = pZ m ' (r + v/2m(/vo - 1) - y/2m(N-l)), 



(G) 
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Fig. 1. — (a) Single-particle densities for 15, 20 and 25 particles in the v = 1/3 quantum Hall regime, 
(b) The same densities as in (a), but relatively shifted (see text). 



one finds, after a little algebra, that for short-range potentials and large TV the fluctuations 
in A2 scale like < (Ag ) 2 >~ N~ 3/2 . Here < ... > denotes average over realizations of the 
random potential. In particular, for delta- function potentials of density no and typical strength 
Vq, one finds < (A^ N) ) 2 >= A m n V^ TV -3 / 2 , where A m is TV-independent for large TV. In fig. 2a 
we show the averaged PSF, < (A^ ) 2 >, in the v = 1, 1/3 and 1/5 quantum Hall regimes. In 
agreement with the above argument, the average PSF scale with TV~ 3/ " 2 (continuous curve). In 
fig. 2b we plot the numerically calculated A m (scaled by the large- TV value A\ — > l/87r 3/2 ) as a 
function of TV. As expected, it is indeed TV-independent within the accuracy of the calculation. 
The numerically deduced values lead to 

< (A< N) ) 2 > 1/3 / < (A< N) ) 2 >i = 5.0 ±0.3 

< (A< N) ) 2 > 1/5 / < (A<">) 2 >i = 9.2 ±0.5. (7) 

These calculations clearly demonstrate the enhancement of the PSF due to the increased role 
of correlations, as the filling factor is lowered [fi0[ . (Note that the PSF are independent of the 
charging energy in this regime). 

As the range of the random potentials, a, increases, the PSF are reduced and the ratios 



4 




10 15 
Number of particles 



Fig. 2. - (a) Averaged peak-spacing fluctuations for v = 1, 1/3, and 1/5. The continuous lines are 
fits to A m /N 3 / 2 (y = 1/3 and 1/5) and the exact result for v = 1. (b) The same data as in (a), 
normalized by the large- iV v = 1 exact result. 



(eq.(0)) decrease towards unity (fig. 3). However, once a becomes of the order of the size of 
the dot, R oc V~N, the PSF become independent of system size and the ratio increase again, 
leading to a nonmonotonic dependence on a. In the limit of very smooth potentials, a — > oo, 
one can show that the PSF are given by 

< (AH 2 >= 5 " < ^ > m (/ d\ r 2 A 2 p^(r)) 2 • (8) 

The quantity in parentheses is twice the second derivative of the total angular momentum 
with respect to N and is equal to 2m. Thus in this limit, one finds 

< (AH 2 > 1/3 / < (A<"») 2 >i = 9 

< (A^) 2 >i/ 5 / < (A< N) ) 2 >i = 25, (9) 

namely there is an enhancement of the PSF also in this limit. As the range of the potentials 
is determined experimentally by the distance to the donor layer, which is of the same order 
of the dot size, the ratio a/R can be varied experimentally in the relevant regime and these 
predictions, including the nonmonotonicity, could be tested. 

The observed features can be understood as follows. Based on Wen's description of the 
edges of the quantum Hall liquid pT[, Kinaret et al. have shown that for filling factors 
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Fig. 3. - The ratio of the average peak-spacing fluctuations between the v = 1 and the v — 1/3 
quantum Hall regimes as a function of the range of the potentials, scaled by the dot size R oc y/N. 
inset: The distribution of the peak-spacing fluctuations for two densities of random potentials. For 
the larger densities the distribution approaches a Gaussian. 



v = l/m, A 2 = mv/R, where R is the radius of the system and v the edge velocity. This 
result can be simply understood - when the number of particles increases by one, the energy 
of the highest-occupied angular momentum state increases by WdR/dN, where VV^ is the 
potential gradient near the edge of the system. Since R ~ \/2mN this immediately leads 
to A 2 ~ mW/R oc mv/R. This result implies < (A 2 ) 2 >~ m 2 < (W) 2 > /N. As 
< (Vy) 2 >^ n < Vq > /er 6 , where n is the number of random potentials, one immediately 
reproduces eq.(^). As the derivation here relied on the sharpness of the edge it only strictly 
applies in that limit. For finite-range potentials one expects the internal structure of the edge 
to modify this result, which in indeed observed numerically. 

Lastly we consider the full distribution of the PSF. For simplicity, we concentrate on 
delta-function potentials. In this case the distribution is defined by 

P(A 2 ) = j Y[d 2 n dVi P(ri,Vi) s |a 2 - VV^pCnH , (10) 

where P{ri,Vi) is the probability of finding a delta-potential of magnitude Vj at point r^. 
Assuming a uniform distribution of potentials of density no, and a Gaussian distribution of 
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their amplitudes, with width Vo, one finds 

p (^) = / II e ^ /(2E2) ; s - v o E ( A ^)) 2 • (ii) 

When the potentials are dense, the sum (in the expression for E) can be approximated by 
an integral, and we immediately find that the A2 is distributed normally with a width that 
has been calculated above. Note that this derivation did not involve any information about 
the quantum Hall system (except the assumption that one can treat the potential fluctuations 
perturbatively), and is just a manifestation of the law of big numbers. 

When the density of the potentials is lowered, there is a finite probability that none of the 
potentials affects the region of interest - the edge of the system. This enhances the probability 
that A2 = 0. In the inset of fig. 3 we plot the distribution of the PSF for two different densities 
of potentials. We indeed see that for low enough density there is a sharp peak at A2 = 0, 
while the distribution approaches a Gaussian for large enough density. We expect a similar 
effect when the range of the potential increases. 

To conclude, we have demonstrated that the increased importance of the correlations in 
the quantum Hall regime leads to an enhancement of the peak-spacing fluctuations. This can 
be attributed to the increased rigidity of the ground-state wavefunction as the filling factor 
becomes lower. 
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